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Evaluating Thermal Fire Detection Systems [English Units] 


David W. Stroup 
David D. Evans 
Phyllis Martin 


ABSTRACT 


This report presents a methodology for evaluating heat detection systems 
installed in buildings. Previous work for use primarily in designing new 
thermal fire detection systems was used as a starting point. The previous work 
was enhanced and supplemented to make it more useful for evaluating existing 
systems. The resulting equations were programmed into a user interactive 
computer program. This program is available in both BASIC and FORTRAN and will 
run on mainframes as well as personal computers. In addition, a modified 
version of the FORTRAN program was used to develop an extensive set of tables 
listing detector activation times for given building geometries, detector 
characteristics, and fire growth rates. These tables are useful for quick 
evaluation of alternative heat detector installations. Finally, practical 
examples are included to illustrate the use of the tables and computer 
programs. 


Keywords: fire alarm systems; fire detection; fire detection systems; fire 
hazard assessment; fire protection; fire suppression; heat detectors; sprinkler 
systems. 


1. INTRODUCTION 


Studies of the response of heat detectors to fire driven flows under 
confined ceilings have been conducted since the early 1970's [in 2 aoa sito 
Results of these largely experimental studies have been used to develop 
correlations of data that are useful under a broad range of fire conditions and 
building geometries. These correlations have been used to construct 
engineering methods to determine heat detector spacing, sprinkler response 
time, and smoke detector alarm times for industrial buildings where large 
undivided ceilings over storage and manufacturing facilities are common. The 
National Fire Protection Association (NFPA) has adopted an alternate design 
method for determining heat detector spacing which is based on some of these 


correlations. 


This alternate design method is contained in Appendix C of the NFPA 72E 


' which consists of 


standard [6], "Guide for Automatic Fire Detector Spacing,' 
tables listing maximum allowable detector spacings for given room sizes and 
fire growth rates calculated from correlations of experimental measurements. 
The tabular form of presentation is useful for specifying design requirements 
in new installations. However, the tables are inconvenient to use in 
evaluating changes to existing detector systems where the spacing is fixed. 
There appears to be a need for a method which determines detector activation 


times for specified building geometries, detector characteristics, detector 


spacings, and fire growth rates. This method would be useful in evaluating the 


1 Numbers in brackets refer to literature references at the end of this repor 


Z 


expected fire detection time for existing and proposed systems for use in 
determining available safe egress time, and the relative effectiveness of 


different detectors used at the same spacing. 


This paper presents the results of work to develop an alternative method 
appropriate to existing installations. As part of this study, the basis for 
the calculation method published in Appendix C of the NFPA 72E standard was 
determined. Alternate correlations of the experimental data, used to develop 
the results presented in Appendix C of NFPA 72E, have been used to construct a 
FORTRAN program (DETACT-T2). This computer program evaluates the response time 
of heat detector systems [7]. In this work, a modified version of the 
DETACT-T2 program was used to generate tables of detector response times for a 
wide range of building geometries, detector characteristics, and fire growth 
rates. These calculated values for response time agree to within 5 percent of 


those published in the tables contained in Appendix C of the NFPA 72E standard. 


2. DESCRIPTION OF THE MODEL 


The problem of heat detector response is solved by predicting the time 
dependent temperature of the detector sensing element up to the point when it 
is heated to the specified alarm conditions. In order to predict the 
temperature of a heat detector element, it is necessary to describe the fire 
and the fire generated environment to which the detector is exposed. The 
primary fire related parameters which influence the response of heat detectors 


are the gas temperature and velocity. Once the fire-generated environmental 


characteristics are known, a model of the response of the thermal sensing 
element in the detector is employed to calculate the detector element 


temperature as a function of time. 


2.1 Characteristic Fire 


Extensive full-scale tests have been conducted to determine the response 
of heat detectors [1-4]. These tests have studied two classes of fires. The 
first of these are steady fires. These fires have energy release rate 
histories which do not change with time. Results from work using this type of 
fire may be used to calculate heat detector response if the quasi-steady state 
assumption is invoked. In this case, a varying energy release rate history is 
approximated as a series of steady fires. In this type of analysis, the 
transport time of the gases to the detector is ignored. This usually results 
in a predicted actuation time which is earlier than that encountered in actual 


practice. This may or may not be acceptable depending on the case. 


Another class of fires for which research results are available are called 
wesequaneay (t?) fires. These fires have energy release rate histories which 
increase proportionally with the second power of time from ignition. This 
class of fires was used to develop the tables in the NFPA 72E document and will 


be used in this paper. 


Appendix C of the NFPA 72E standard lists three fires which have t? type 


fire growth behavior. A slow developing fire is defined as one which would 


take 600 seconds (10 minutes) from the time of flaming ignition until it 
reaches a heat release rate of 1000 BTU/s (1055 kW). A medium growth rate fire 
reaches a heat release rate of 1000 BTU/s (1055 kW) approximately 300 seconds 
(5 minutes) after flaming ignition occurs. The fast developing fire requires 
150 seconds (2.5 minutes) to reach a heat release rate of 1000 BTU/s (1055 

kW). In addition, a fourth t? fire has been developed in this report. This 
fire has been defined as one which takes 75 seconds (1.25 minutes) to reach a 
heat release rate of 1000 BTU/s (1055 kW). These four fire growth curves are 


plotted in Figure 1. 


Table 1 is a reproduction of a table in NFPA 72E. It lists various 
warehouse commodities and their maximum fire heat release rates per unit floor 
area covered by the commodity. The time required for a fire, occurring in a 
particular commodity, to reach a heat release rate of 1000 BTU/s is also 
listed. If possible, each commodity is identified with a specific fire growth 
behavior - slow, medium, fast, or ultrafast fire. Table 2 is a reproduction of 
another table found in Appendix C of NFPA 72E. It lists the maximum heat 
release rates for various furnishing type fuel items. Both of these tables are 
useful for estimating fire growth rates that may be expected from various fuel 


packages. 


2.2 Temperature and Velocity Correlations 


Values of the time dependent gas velocity and temperature for t? fires may 


be obtained from work by Heskestad and Delichatsios [4]. From experiments and 


use of dimensional arguments, they developed correlations of maximum 
temperature rise and velocity for unconfined ceiling layer flows as a function 
of radial distance from the axis of fire plume impingement. These correlations 


of dimensionless maximum temperature rise and dimensionless maximum velocity 


are: 
ou %* 
0, ty < (to)¢ 
aT, = (1) 
% 4/3 
ey = O505) Cho) oa 
ee %* % 
0.188 + 0.313 r/H >» ac (t5)¢ 
Unf JAD, = 0659 (city (2) 
where 
< Sep eS us ls 
ue Sr ee te es 
res ae nin  GuiaD eel? a=] 
A = G@ © 
go) iC ao 
a = Q/ ee 
* 
(to). = 0.954 (1 + r/H) 


These equations may be used to predict the fire-generated environment (gas 


temperature and velocity) near a detector. 


2.3 Detector Response Model 


Characterization of the thermal response of heat detector and sprinkler 
thermal sensing elements is discussed by Heskestad and Smith [8], and Evans 
[9]. The model for the detector sensing element temperature is based on a 
convective heat transfer process. The first order differential equation that 


describes the rate of temperature increase of the sensing element is: 


aie cea: Bel (3) 


This equation predicts the temperature change of a detector with known thermal 
characteristics to a given time dependent gas flow temperature and velocity. 
In particular, the rate of change of the detector element temperature at any 
time is equal to the gas velocity to the one-half power multiplied by the 
difference between the temperature of the hot fire generated gas flow near the 
detector and the temperature of the detector element. All divided by the 


detector RTI or response time index. 


The RTI is a measure of the detector's thermal lag and can be determined 
for sprinkler or heat detector sensing elements using the plunge test [8]. 
NFPA 72E reports detector sensitivity in terms of time constants, t. At any 
given time, the time constant for a given type of detector may be calculated 


from the RTI using: 


[2 


Te RE of ut (4) 


Tables 3 and 4 relate heat detector listed spacings to RTI values for fixed 
temperature detectors and rate of rise detectors respectively. These Tables 


were derived from results published in references 3 and 5. 


Beyler [10] used the correlations of Heskestad and Delichatsios [4] to 
derive a new form of equation (3). He was able to solve this new form of the 
heat detector response equation analytically. Using the condition that the 
initial detector element temperature is equal to ambient temperature Ge = Le 
the equations for calculating the response of fixed temperature and rate of 


temperature rise detectors are, from Beyler: 





* * 1/4 
dT (4/3) (AT/AT.) (AT.) 
s Ky ee Giese) (5) 
dt (t/t,) (0.188 + 0.313 r/H) 
T = 7.(0) = .CAT/AT") AL. | te (Glee) (6) 
s Ss a 2 2 Y 


where 


¥ ¥ 
U2 U AT 
la ay ies 2 ead tt 
Y = i x i at a (0.188 + 0.313r/H) 
Uy AT. to 


Using equations (1) and (2), equations (5) and (6) may be solved to predict the 
response of sprinkler and heat detector sensing elements to fire-driven ceiling 


jet flows. 


2.4 Velocity Correlation for r/H < 0.3 


While this exact solution to equation (3) is applicable to any r/H value, 
the dimensionless velocity correlation, previously discussed, is only valid for 
r/H values greater than 0.3. In an effort to extend the usefulness of the 
detector response calculations, an analysis was made of the available data for 
ceiling jets at r/H values less than 0.3. Results of this analysis indicated a 
lack of much quantitative information for this region. Therefore, an analysis 
to develop a rational estimation of velocities in the ceiling-jet close to the 
plume axis was undertaken by studying a simple extension of the velocities in 


the plume. 


The dimensionless velocity and temperature change relations proposed by 
Heskestad and Delichatsios [4] together with equations proposed by Zukoski 
et. al. [11] to describe fire plume maximum temperature and velocity 


distributions were taken as a starting point. 3 


1/5 tls yl/5 


U u/ [A ] (7) 
AT, = ot / tat’ (T_/g) CEs e wali (8) 
Mote emt co (GO)? (9) 
u = G gil? yl/2 (g%y1/3 (10) 


These equations were used to obtain an expression of dimensionless maximum gas 


velocity in the same form as equation (2). 


Ua! Jat, = Gey IC (11) 


Using values suggested by Zukoski et. al. [11], the dimensionless gas velocity 
is found to be a constant at r/H values less than 0.3. The complete 


dimensionless gas velocity equation for all values of r/H is: 


i = S87 |) Jems, ell S Oo3 
uU, / Jar, = entes (12) 
O59 Gal) efit > Oo 


This equation is plotted in Figure 2. The two curves merge nicely. This 
equation together with the previously discussed correlations and analytic 
solution may be used to calculate any aspect of the response of thermal 


detectors to t? fires. 


3. IMPLEMENTATION OF THE MODEL: COMPUTER PROGRAMS 


3.1 General Program Structure 


The equations discussed in Section 2 were programmed into a user 
interactive computer code called DETACT-T2. To evaluate the response of a 
detector, the user enters values for ambient air temperature, detector 
activation temperature and rate of temperature rise, detector RTI, fuel to 
ceiling distance, detector spacing, and a fire growth rate constant, a (for t? 
fires). The program solves the appropriate equation to obtain the response 


time for a fixed temperature detector and a rate of rise detector with the 


10 


given characteristics. The fire energy release rate at each time step is also 


calculated. 


The computer program is available in two versions. One version, contained 
in Appendix A, is written in FORTRAN 77. The program contains approximately 
400 lines of source code and requires about 64k of computer memory in compiled 
form for execution. The source code may be compiled on any computer supporting 
an ANSI standard FORTRAN 77 compiler. The program has been compiled and run 
successfully on IBM* and compatible personal computers operating under 
MS-DOS?. A specially modified version of this program was used to generate the 


tables discussed in section 4. 


The other version of the program, listed in Appendix B, is written in 
BASIC. It consists of approximately 350 lines of source code. This program 
has also been successfully tested on IBM and compatible micro-computers. Both 
programs have essentially the same structure and use the same variable names. 


To insure convergence, all variables are in double precision. 


The program consists of three major sections. The main program reads and 
lists the required input data, performs the calculations necessary to set up 
the equations, calls the solution subroutines, and writes out the results. The 
second section is a subroutine which uses a Newton-Raphson technique to solve 
the equation describing the activation time of a fixed temperature detector. 

2 IBM is a trademark of International Business Machines Corporation and 
does not represent an endorsement by the National Bureau of Standards. 


3 MS and Microsoft are trademarks of Microsoft Corporation and do not 
represent an endorsement by the National Bureau of Standards. 


it 


The last section is a subroutine which uses a bisection method to solve the 
rate of temperature rise detector activation time equation. Each of these 


sections is described in more detail below. 


3.2 Main Program 


3.2.1 Data Input Section 


This section of the main program requests information from the user 
concerning the building geometry, the detector characteristics, and the fire 
growth rate. The user may also elect to enter data in either English or 


S.1. units. 


3.2.2 Calculation Initialization 


Several preliminary calculations must be made before the activation time 
equations may be solved. These calculations are done in the second section of 
the main program. In this section, the dimensionless time, velocity, and 
temperature are calculated. Based on the user-specified detector spacing, the 
maximum radial distance of the detector from the fire is calculated. After 
these calculations have been completed, the two activation time calculation 
subroutines are called. The results from these subroutines are converted into 


activation times, and the heat release rates are calculated. 
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3.2.3 Conversion and Output Section 


The last section of the main program writes out the solutions to the 


equations. The results are printed in both English and S.I. units. 


3.3 Activation Time Equation Solution: Fixed Temperature 


The equation describing the activation time of a fixed temperature 
detector cannot be solved explicitly for activation time. Therefore, a 
numerical technique must be used. A Newton-Raphson technique was selected for 
use in solving the fixed temperature detector equation. This technique is 


discussed in most numerical methods textbooks [12]. 


The Newton-Raphson (or simply Newton's) method is one of the most powerful 
and well-known numerical methods for finding a root (solution) of f(x) = 0. To 
use this routine, the activation time equation is first adjusted to fit the 
appropriate form. This method converges very rapidly for the cases considered 


here. 
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3.4 Activation Time Equation Solution: Rate of Temperature Rise 


A bisection technique was selected for use in solving the rate of rise 
activation time equation. This technique is also discussed in most numerical 
methods textbooks [12]. It was chosen for use in this instance because the 
Newton method would have introduced stability problems into the calculation. 


This method also converges very rapidly for the cases considered here. 


4. IMPLEMENTATION OF THE MODEL: TABLES OF RESULTS 


The computer programs calculate detector activation times one case at a 
time. In many instances, it is necessary to evaluate a number of 
alternatives. It would be tedious and inconvenient to run a computer program 
to examine each case of interest. In light of this fact, the FORTRAN version 
of the computer program was used to generate an extensive set of tables. These 
tables list detector activation time as a function of building geometry, 


detector characteristics, and fire growth rate constants. 


Each entry in the tables was checked by first calculating the response 
time of the given detector numerically. Then the calculated time is 
substituted into the exact form of the equation, and it is solved for the 
temperature of the detector sensing element at that time. In each case, the 
calculated temperature agreed with the originally assumed one to within an 


accuracy of 0.005 percent. 
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The tables cover a wide range of parameter variations. Tables have been 
developed for each of the four fire growth rates that were described earlier. 
Tables are available in both English and S.I. units. The English units tables 
are contained in Appendixes C and D of this report. The S. I. units tables are 


available in the metric version [13]. 


For use with English units, tables cover detector spacings of 0 to 60 feet 
and fuel to ceiling heights of 4 to 80 feet. RTI values range from 0 to 800 
(ft-s)!/?2 for fixed temperature detectors and from 100 to 2000 (ft-s)!/? for 
rate of temperature rise detectors. Temperature differences for fixed 
temperature detectors range from 10 to 140 degrees F, and available rate of 


temperature rise values are 15, 20, and 25 degrees F/minute. 


The tables of detector activation times are organized according to: type 
of detector (fixed temperature or rate of temperature rise), detector spacing, 
and type of fire (ultrafast, fast, medium, or slow). Appendix C contains the 
tables listing activation times for fixed temperature heat detectors, and 
Appendix D contains the tables listing activation times for rate of temperature 
rise detectors. In each appendix, the tables are arranged in order of 
increasing detector spacing (i.e. 0 ft, 2 ft, 4 ft, etc.). For each spacing, 
the tables are arranged according to fire type. Ultrafast is first, followed 


by fast, medium, and slow. 


The asterisks next to some of the table entries indicate that the detector 
location places it within the r/H < 0.3 region. This region is where the 


theoretical dimensionless velocity correlation, discussed in section 2.4, is 
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applicable. These values are flagged because the calculation of gas velocity 
in this region is based primarily on a theoretical analysis and not on a 
combination of theory and experiment. In addition, this region is not 


considered in the NFPA 72E Appendix C. 


Finally, just as in Appendix C of NFPA 72E, for the tables in this report, 
the fire plume is assumed to impinge on the ceiling at the center of the square 
array of detectors spaced as indicated. So for detectors spaced at 10 
ft (3.048 m), the distance of each detector from the gas impingement point 
would be 7.07 ft (10 /2/2). Each table indicates this distance in parenthesis 


after the table spacing value. 


5. EXAMPLES 


5.1 Using the Tables 


5.1.1 Example 1: Comparison with NFPA 72E Table Results 


The first example will show that these new tables agree quite well with 
the results from the tables in Appendix C of NFPA 72E even though each set of 
tables were generated using different correlations of experimental data. 

Assume a system of detectors is installed using a 14 foot spacing in a room 
with a 12 foot ceiling height. A medium fire growth rate is expected. The 
detectors have an RTI rating of 50 (ft-s)1+/? which corresponds to a 22.4 second 


time constant. The ambient temperature is 70 degrees F and the rated 
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activation temperature of the detector is 190 degrees F. This gives a 
temperature difference of 120 degrees. 


AT = Tact - Tamb = 190 - 70 = 120 degrees F 


Looking in Appendix C, a detection system installed using a 14 foot 
spacing would detect a medium growth rate fire 4.91 minutes after ignition. If 
we next look for similar conditions in the 72E tables [NFPA, p. 72E-36], we 
find that a detector installed under these conditions with a time constant of 
25s will activate when a medium growth rate fire reaches a heat release rate of 
1000 BTU/s. A medium growth rate fire will reach this heat release rate 300 


seconds or five minutes after ignition. 


The detector listed in the NFPA table has a time constant of 25s. This is 
slightly higher than the 22.4s for the detector listed in the NBS table. 
Therefore, it is expected that the detector listed in the NBS table would have 
an activation time slightly lower than the one in the NFPA table. After 
extensive cross checks between the new NBS tables and the 72E tables, the 


results were found to agree to within 5 percent. 


5.1.2 Example 2: Fixed Temperature vs. Rate of Rise 


This example will show how alternative heat detector installations may be 
compared using the tabular results. In particular, the differences between 
fixed temperature and rate of rise detectors will be compared. Suppose a 


detection system is installed in a room with a ceiling height of 40 feet and 
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detection of a medium fire is desired. The detectors are spaced 20 feet 

apart. Currently, fixed temperature detectors with an RTI of 200 (ft-s)}/2 and 
an activation temperature rating of 170 degrees F are installed in the room. 
Ambient temperature is 70 degrees F. (Temperature differential of 100 
degrees.) Will the use of rate of rise detectors with the same RTI as the 
fixed temperature detectors and an activation rate of rise rating of 15 degrees 
F/min yield a lesser time to activate a detector in response to a medium growth 


rate fire? 


From Appendix C, the activation time of the fixed temperature detectors, 
installed using a 20.0 foot spacing and exposed to a medium growth rate fire, 
is 12.27 minutes. The time required to activate the rate of rise detectors is 
obtained from Appendix D. The activation time for these detectors is 15.29 
minutes. From Figure 1, the heat release rate at activation of the rate of 
rise detectors would be 9342 BTU/s compared to 6016 BTU/s for the fixed 


temperature detectors. 


It would take about three minutes longer to activate the rate of rise 
detectors, and the fire heat release rate would have increased by about 3326 
BTU/s. This is interesting in light of the generally held belief that rate of 
rise detectors activate before fixed temperature detectors. This is not true 
when the detectors are located in areas with high ceilings and the fire does 
not grow rapidly. This example illustrates how quick and sometimes surprising 
comparisons may be made using the tabular form of the detector activation time 


results. 
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5.2 Example 3: Using the Computer Program 


The last example will be used to show how the computer programs work. 
Both the BASIC version and the FORTRAN version ask similar questions and 
operate in a similar manner so this discussion is applicable to both programs. 
The problem will be to calculate the activation time for fixed temperature and 
rate of rise detectors installed, using a 10 foot spacing, in an area with a 
ceiling height of 12 feet. The detectors have an RTI of 670.8 (ft-s)1/*. The 
detector activation temperature is 130 degrees F, and the activation rate of 
rise is 15 degrees F/min. Ambient temperature is 70 degrees F. (User input is 


preceded by a > and, where possible, typed in lower case.) 


ENTER 1 FOR ENGLISH UNIT INPUT 
2 FOR METRIC UNIT INPUT 
>1 


ENTER THE AMBIENT TEMPERATURE IN DEGREES F. 
>70 


ENTER THE DETECTOR RESPONSE TIME INDEX (RTI) IN (FT-SEC)**1/2. 
>670.8 


ENTER THE DETECTOR ACTIVATION TEMPERATURE IN DEGREES F. 
>130 


ENTER THE DETECTOR RATE OF RISE IN DEGREES F/MIN. 
lS) 


ENTER THE CEILING HEIGHT IN FEET. 
N72 


ENTER THE DETECTOR SPACING IN FEET. 
>10 


ENTER: S FOR SLOW FIRE GROWTH RATE, 
M FOR MEDIUM FIRE GROWTH RATE, 
F FOR FAST FIRE GROWTH RATE, 
U FOR ULTRAFAST FIRE GROWTH RATE, OR 
O FOR OTHER. 
>m 
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RESULTS: 


CEILING HEIGHT 


= 12.0 FEET ( 3.66 METERS) 
DETECTOR SPACING = 


10.0 FEET ( 3.05 METERS) 
DETECTOR RTI = 670.8 (FT-SEC)**1/2 ( 370.3 (M-SEC)**1/2) 
FIRE GROWTH RATE CONSTANT = .1111E-001 BTU/(SEC**3) 
(.1171E+002 JOULES/SEC**3) 
FOR TEMPERATURE ACTUATED DETECTOR: 
ACTIVATION TEMPERATURE = 130.0 DEGREES F (54.4 DEGREES C) 
TIME OF ACTIVATION = 4.96 MINUTES 


HEAT RELEASE RATE = .9859E+003 BTU/SEC 
(.1040E+004 KILOJOULES/SEC) 


FOR RATE OF RISE ACTUATED DETECTOR: 


ACTIVATION RATE OF RISE = 15.00 DEGREES F/MIN 
( 8.33 DEGREES C/MIN) 


TIME OF ACTIVATION = 3.04 MINUTES 


HEAT RELEASE RATE = .3707E+003 BTU/SEC 
(.3911E+003 KILOJOULES/SEC) 


The results show that the fixed temperature heat detector would activate 
approximately 5 minutes after the fire reaches a flaming state. The heat 
release rate at this time would be 986 BTU/s. A rate of rise detector would 
activate at about 3 minutes with a corresponding heat release rate of 371 
BTU/s. (If metric units had been selected, the input requests would have 


called for data in metric units instead of English units.) 
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6. SUMMARY 


This report presented a methodology for evaluating heat detection systems 
installed in existing buildings. Previous work for use primarily in designing 
new systems was enhanced and supplemented to broaden its application. The 
resulting equations were programmed into a user interactive computer program. 
This program is available in both BASIC and FORTRAN and will run on mainframes 
as well as personal computers. In addition, a modified version of the FORTRAN 
program was used to develop an extensive set of tables listing detector 
activation times for given building geometries, detector characteristics, and 
fire growth rates. These tables are useful for quick evaluation of alternative 
heat detector installations among other uses. Practical examples were utilized 


to illustrate the use of the tables and computer programs. 
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. Wood pallets, stack 


Settee hi (6-12% 
moisture) 


- Wood pallets, stack 


5 ft high (6-12% 
moisture) 


- Wood pallets, stack 


10 ft high (6-12Z% 
moisture) 


- Wood pallets, stack 


16 ft high (6-122 
moisture) 


. Mail bags, filled, 


stored 5 ft high 


. Cartons, compart- 


mented, stacked 15 
ft high 


Paper, vertical rol 
stacked 20 ft high 


. Cotton (also PE, 


PE/Cot, Acrylic/ 
Nylon/PE), gar- 
ments in 12 ft high 
rack 


. Cartons on pallets 


rack storage, 15-30 
ft high 


Adapted from Table 


TABLE 1. WAREHOUSE MATERIALS* 


Maximum 
Heat 
Time to Release 
Reach per Unit 
1000 BTU/s Floor Area 
(sec) (BTU/s/ ft?) 
150-310 110 
90-190 330 
80-110 600 
75-105 900 
190 35 
60 200 
ls, 
15-28 = 
20-42 - 
40-280 - 


C-2-2.2.1(a) of Ref. 5 
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Classification 
(s - slow) 
(m - medium) 
(£ - fast) 

(u - ultrafast) 


m-f 


10. 


Mate 


12. 


3G 


14. 


NS)o 


16. 


IZ) 


18. 


Qe, 


20; 


TABLE 1. WAREHOUSE MATERIALS (continued 


Maximum 
Heat 
Time to Release 
Reach per Unit 
1000 BTU/s Filcor Area 
(sec) (BTU/s/ft?) 


Paper products, 

densely packed in 

cartons, rack 

storage, 20 ft high 470 = 


PE letter trays, 
filled, stacked 5 ft 
high on cart 190 750 


PE trash barrels in 
cartons stacked 15 ft 
high 55 250 


FRP shower stalls in 
cartons, stacked 15 
ft high 85 110 


PE bottles packed in 
Item 6 85 550 


PE bottles in car- 
tons, stacked 15 ft 
high 75 170 


PE pallets, stack 3 ft 
high 130 = 


PE pallets, stack 6-8 
ft high 30555 = 


PU mattress, single, 
horizontal 110 - 


PF insulation board, 
rigid foam, stacked 
ISeftehigh 8 170 


PS jars packed in 
Item 6 55 1200 
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Classification 
(s - slow) 
(m - medium) 
(f£ - fast) 

(u - ultrafast) 


Zi 


22. 


23. 


24. 


25% 


26. 


27. 


28. 


29. 


30. 


31. 


TABLE 1. WAREHOUSE MATERIALS (continued 


PS tubs nested in 
cartons, stacked 14 
ft high 


PS toy parts in car- 
tons, stacked 15 ft 
high 


PS insulation board, 
rigid, stacked 14 ft 
high 


PVC bottles packed 
in Item 6 


PP tubs packed in 
Item 6 


PP and PE film in 
rolls, stacked 14 ft 
high 


Distilled spirits in 


barrels, stacked 20 ft 


high 

Methyl alcohol 
Gasoline 
Kerosine 


Diesel Oil 


Abbreviations: 


Time to 
Reach 
1000 BTU/s 
(sec) 


105 


110 


10 


40 


FRP - Fire Retardant Plastic 


= Polyethylene 


- Phenolic Resins 


a Polypropylene 


25 


Maximum 
Heat 
Release 
per Unit 
Floor Area 
(BTU/s/f£t?) 


450 


180 


290 


300 


390 


350 


Classification 
(s - slow) 
(m - medium) 
(£ - fast) 

(u - ultrafast) 


Polysytrene 
Polyurethane 
Polyvinyl Chloride 
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(ey) 


‘ Adapted 


TABLE 2. MAXIMUM HEAT RELEASE RATES FROM 


FIRE DETECTION INSTITUTE ANALYSIS ~ 


. Medium wastebasket with milk cartons 
. Large barrel with milk cartons 
. Upholstered chair with polyurethane foam 


. Latex foam mattress (heat at room door) 


Furnished living room (heat at open door) 


from Table C-2-2.2.1(b) of Ref. 3 
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Approximate 
Values 
BIU/s 
100 
140 
350 
1200 


4000 - 8000 


TABLE 3. RTI VALUES FOR ANY LISTED DETECTOR’ 
(FIXED TEMPERATURE TYPE) [(ft-s)!/?] 








Listed UL Listed Activation Temperature ne 
Spacing Listed 
(it) 07! 128°F  — 1358RPSeD SOR iw TooeRderami7oerben |196eN lien 
10 894 738 586 436 358 217 436 
15 559 425 349 246 199 101 246 
20 369 302 235 157 116 38 157 
25 277 224 174 107 72 107 
30 BNP 179 136 81 49 81 
40 159 128 92 40 
50 132 98 67 
70 81 54 20 
NOTE: These RTIs are based on an analysis of the Underwriters 


Laboratories and Factory Mutual listing test procedures. Plunge test (see 
Ref. 7) results performed on the detector to be used will give a more 
accurate response time index. 


Abbreviations: 
il Factory Mutual 
UL = Underwriters Laboratories Inc. 


* adapted from Table C-3-2.1.1 of Ref. 5 
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TABLE 4. RTI VALUES FOR ANY LISTED DETECTOR™ 
(RATE OF RISE TYPE ft=s)?72] 


Listed UL Listed Activation Rate of Temp. Rise 
Spacing 
(ft) 15 °F/min 20 °F/min 25 °F/min 
10 1834 1308 984 
7368) 1453 1073 805 
15 1185 872 637 
20 872 581 425 
30 559 380 280 
40 447 291 206 
50 425 246 161 
Abbreviation: 
Wi Underwriters Laboratories Inc. 


* adapted from Ref. 3 
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APPENDIX A. DETACT-T2 
FORTRAN VERSION 
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Qa 


Ailey ar@sari@ai@s@r@rai@rG @iGr@i@t@T@r@l@i@ @Gr@r@r@ area @i@l@r@y@ri@i@r@r@ 


PROGRAM DETACT 


THIS IS PROGRAM DETACT-T2. 

IT CALCULATES THE RESPONSE OF THERMAL DETECTORS, BOTH FIXED TEMPERATURE 
AND RATE OF TEMPERATURE RISE, TO FIRES WHOSE HEAT RELEASE RATE GROWS 
PROPORTIONALLY WITH THE SQUARE OF TIME FROM IGNITION. 


WRITTEN BY DAVID W. STROUP, CENTER FOR FIRE RESEARCH, NBS 


THIS PROGRAM IS A CONTRIBUTION OF THE NATIONAL BUREAU OF STANDARDS AND 
IS NOT SUBJECT TO COPYRIGHT. 


VERSION 1.1 

RHA RRRRRRARRERE NOTATION SHARMA 
INPUT: 

J - UNITS CODE (1 OR 2) 


1 - INPUT DATA IN ENGLISH UNITS 
2 - INPUT DATA IN METRIC UNITS 
TAMB - AMBIENT TEMPERATURE 


RTI - DETECTOR RESPONSE TIME INDEX 

TACT - DETECTOR ACTIVATION TEMPERATURE 

ROR - DETECTOR ACTIVATION RATE OF RISE 

HF - CEILING HEIGHT 

ZF - DETECTOR SPACING 

M - GROWTH FACTOR CODE, CHARACTER VARIABLE (S, M, F, U, OR 0) 
ALPHA - FIRE GROWTH RATE CONSTANT 

OUTPUT: 

at - TIME OF ACTIVATION FOR FIXED TEMPERATURE DETECTOR 
QD - HEAT RELEASE RATE AT TIME OF ACTIVATION, T 

TR - TIME OF ACTIVATION FOR RATE OF RISE DETECTOR 

QDR - HEAT RELEASE RATE AT TIME OF ACTIVATION, TR 


IERR - ERROR CODE (O OR 1) 
0 - SUCCESSFUL 
1 - UNSUCCESSFUL 


9 ale alo ale w%o ale ale ale ate ale ale als ale ales ale alo ale ale ale ale ale alo ale ale ale ale al, Yo we 0! alo 9p ale lo ale alo ale ale als alo ale alo alo ale ale alos lo ale als ale ale alo we als ale ale ale alle als ale ale alo ate Xp ale als alo alo 0! 
Sede se ce Sees ESSER SSRI EE IIE EEE IE EERE EER ERE RERREERER 


IMPLICIT DOUBLE PRECISION (A-H,0-Z) 
CHARACTER M 

DATA IRTTY, IWTTY/5,6/ 

DATA GE,CPE, RHOE/32.2D0,0.24D0,0.0735D0/ 
DATA GM,CPM,RHOM/9.8D0,1.0035D0,1.1768D0/ 
TERR = 0 

WRITE (IWTTY,1) 

FORMAT (' DETACT-T2 VERSION 1.0'/ 

& ' WRITTEN BY D.W. STROUP (1985)'// 
& ' CONTRIBUTION OF THE '/ 
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NATIONAL BUREAU OF STANDARDS (U.S.).'/ 
NOT SUBJECT TO COPYRIGHT. '// 
CALCULATES DETECTOR ACTUATION TIME'/ 
BELOW UNCONFINED CEILINGS WITH '/ 

TIME SQUARED FIRE GROWTH RATES.'//) 


— 2 Bm 


ENTER INPUT DATA 


WaAAN 


CONTINUE 
WRITE (IWTTY,10) 
10 FORMAT (' ENTER: 1 FOR ENGLISH UNIT INPUT'/8X,'2 FOR METRIC UNIT I 
&NPUT' / ) 
READ (IRTTY,*) J 
IF (J.EQ.1) THEN 
WRITE (IWTTY, 20) 
20 FORMAT (' ENTER THE AMBIENT TEMPERATURE IN DEGREES F.'/) 
READ (IRTTY,*) TAMB 
WRITE (IWTTY, 30) 
30 FORMAT (' ENTER THE DETECTOR RESPONSE TIME INDEX (RTI) IN (FT-SEC) 
&**1/2.'/) 
READ (IRTTY,*) RTI 
IF (RTI.LE.0.000001D0) RTI = 0.000001D0 
RTIR = RTI 
WRITE (IWTTY, 40) 
40 FORMAT (' ENTER THE DETECTOR ACTIVATION TEMPERATURE IN DEGREES F.' 
&/) 
READ (IRTTY,*) TACT 
WRITE (IWTTY,50) 
50 FORMAT (' ENTER A DETECTOR RATE OF RISE IN DEGREES F/MINUTE. '/) 
READ (IRTTY,*) ROR 
WRITE (IWTTY,60) 
60 FORMAT (' ENTER THE CEILING HEIGHT IN FEET.'/) 
READ (IRTTY,*) HF 
WRITE (IWTTY, 70) 
70 FORMAT (' ENTER THE DETECTOR SPACING IN FEET.'/) 
READ (IRTTY,*) ZF 
WRITE (IWTTY, 80) 
80 FORMAT c ENTER: S FOR SLOW FIRE GROWTH RATE '/ 
M FOR MEDIUM FIRE GROWTH RATE'/ 
' F FOR FAST FIRE GROWTH RATE'/ 
U FOR ULTRAFAST FIRE GROWTH RATE OR'/ 
' O FOR OTHER'/) 
READ (IRTTY,410) M 


— QR & 


IF ((M.EQ.'s').OR.(M.EQ.'S')) ALPHA = 0.00277778D0 
IF ((M.EQ.'m').OR.(M.EQ.'M')) ALPHA = 0.01111111D0 
IF ((M.EQ.'f').OR.(M.EQ.'F')) ALPHA = 0.0444445D0 
IF ((M.EQ.'u').OR.(M.EQ.'U')) ALPHA = 0.1777778D0 
IF ((M.EQ.'o').OR.(M.EQ.'O')) WRITE (IWTTY,90) 


IF ((M.EQ.'o').OR.(M.EQ.'O')) READ (IRTTY,*) ALPHA 
90 FORMAT (' ENTER THE FIRE GROWTH RATE CONSTANT (ALPHA) IN BTU/SEC/S 
&EC/SEC.'/) 
100 CONTINUE 
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ELSE 
C 
C METRIC UNIT INPUT 
C 
WRITE (IWTTY,110) 
110 FORMAT (' ENTER THE AMBIENT TEMPERATURE IN DEGREES C.'/) 
READ (IRTTY,*) TAMB 
WRITE (IWTTY, 120) 
120 FORMAT (' ENTER THE DETECTOR RESPONSE TIME INDEX (RTI) IN (M-SEC)* 
&*1/2.'/) 
READ (IRTTY,*) RTI 
IF (RTI.LE.0.000001D0) RTI = 0.000001D0 
RTIR = RTI 
WRITE (IWTTY, 130) 
130 FORMAT (' ENTER THE DETECTOR ACTIVATION TEMPERATURE IN DEGREES C.' 
&/) 
READ (IRTTY,*) TACT 
WRITE (IWTTY,140) 
140 FORMAT (' ENTER A DETECTOR RATE OF RISE IN DEGREES C/MINUTE.'/) 
READ (IRTTY,*) ROR 
WRITE (IWTTY, 150) 
150 FORMAT (' ENTER THE CEILING HEIGHT IN METERS. '/) 
READ (IRTTY,*) HF 
WRITE (IWTTY, 160) 
160 FORMAT (' ENTER THE DETECTOR SPACING IN METERS.'/) 
READ (IRTTY,*) ZF 
WRITE (IWTTY, 80) 
READ (IRTTY,410) M 
IF ((M.EQ.'s').OR.(M.EQ.'S')) ALPHA 
IF ((M.EQ.'m').OR.(M.EQ.'M')) ALPHA 
IF ((M.EQ.'£').OR.(M.EQ.'F')) ALPHA = 46.88888889D0 
IF ((M.EQ.'u').OR.(M.EQ.'U')) ALPHA = 187.5555556D0 
IF ((M.EQ.'o').OR.(M.EQ.'O')) WRITE (IWTTY,170) 
IF ((M.EQ.'o').OR.(M.EQ.'O')) READ (IRTTY,*) ALPHA 
170 FORMAT (' ENTER THE FIRE GROWTH RATE CONSTANT (ALPHA) IN WATTS/SEC/ 
&SEC.'/) 
180 CONTINUE 
ENDIF 


2.930555556D0 
11.72222222D0 


C 
GC wvrererererederererevereresedereveres CAT CULATIONS wevererere vere re voce ve ve ve ve ve ve nk seve ve 26 He Be Fe Be B0 He Bee He He He He He He Fe 
C 


R = 0.5D0*DSQRT(2.0D0)*ZF 
ROH = R/HF 

ROR = ROR/60.D0 

IF (J.EQ.1) THEN 


TAMB = TAMB + 460.D0 
TACT = TACT + 460.D0 
A = GE/(CPE*TAMB*RHOE ) 
G = GE 


TAMB = TAMB + 273.D0 
= TACT + 273.D0 
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C 


C kx PRINT OUT RESULTS “ecececsdede deve ce ve te ve ve ve ve te seve ve se ve de se se se seve 


C 


190 


C 


A = GM/(CPM*TAMB*RHOM*1000.D0) 
G = GM 
ENDIF 


TOTS2 = A**(-1.D0/5.D0)*ALPHA**(-1.D0/5.D0)*HF**(4.D0/5.D0) 
DLTODL = A**(2.D0/5.D0)*(TAMB/G)*ALPHA**(2.D0/5.D0) 


& *HF**(-3.D0/5.D0) 
UOUS2 = A**(1.D0/5.D0)*ALPHA**(1.D0/5.D0)*HF**(1.D0/5.D0) 
DELTD = TACT - TAMB 


IF (ROH.GT.0.3D0) THEN 

UODLTH = 0.59D0*ROH**(-0.63D0) 

ELSE 

UODLTH = 3.87D0/(9.115D0**0.5D0) 

ENDIF 

TS2F = 0.954D0*(1.D0+ROH) 

A2 = (4.D0/3.D0)*DLTODL*UOUS2**(-0.5D0)*UODLTH**(-0.5D0)*RTI/ 
& (TOTS2*(0.188D0+0.313D0*ROH) ) 

C = 1.D0+DELTD/A2 

CALL NWIN (C,Y,IERR, IWTTY) 

IF (IERR.EQ.1) GO TO 380 

DELTS2 = (4.D0/3.D0)*U0US2**(-0.5D0)*UODLTH**(-0.5D0)*RTI*Y/ 
& (TOTS2*(0.188D0+0.313D0*ROH) ) 

TS2 
TS2 = TS2 + TS2F 

T = TS2*A**(-1.D0/5.D0)*ALPHA**(-1.D0/5.D0)*HF**(4.D0/5.D0) 
QD = ALPHA*T**2 

IF (RTIR.LE.2.0D0) RTIR = 2.0D0 
Dl = (4.D0/3.D0)*DLTODL/ (TOTS2*(0.188D0+0.313D0*ROH) ) 
D2 = (3.D0/4.D0)*UOUS2**0 . 5DO*UODLTH**0 . 5D0*(1.DO0/RTIR)*TOTS2 
& *(0.188D0+0.313D0*ROH) 

CALL BISECT (D1,D2,ROR,DELTS2, TERR, IWTTY) 

IF (IERR.EQ.1) GO TO 380 

TS2R = 0.954D0*(1.D0+ROH)+(0.188D0+0. 313D0*ROH)*DELTS2 
& **(3.D0/4.D0) 

TS2R = TS2R + TS2F 

TR = TS2R*A**(-1.D0/5.D0)*ALPHA**(-1.D0/5.D0)*HF**(4.D0/5.D0) 
QDR = ALPHA*TR**2 


* 
* 
a 


ROR = ROR*60.D0 

WRITE (IWTTY, 190) 
FORMAT ('1 RESULTS:') 
IF (J.EQ.1) THEN 


C ENGLISH (METRIC) OUTPUT 


C 


200 


HF2 = HF * 0.3048D0 
ZF2 = ZF * 0.3048D0 
WRITE (IWTTY,200) HF,HF2,ZF,ZF2 


FORMAT ('0 CETLING HEIGHT " = ' 6.2," FEET (',F6.2)!°METERS) '/ 
= (CEO. 20 REET (USRG. 2, )AMEMERS )' ) 


& 5X, 'DETECTOR SPACING 
RTI2 = RTI * (0.3048D0**0.5D0) 


ey 


0.954D0*(1.DO+ROH)+(0.188D0+0. 313D0*ROH)*DELTS2**(3.D0/4.D0) 


WRITE (IWTTY,210) RTI,RTI2 
210 FORMAT ('0 DETECTOR RTI = ',F6.1,' (FT-SEC)**1/2 (',F6.1,' (M-S 
&EC)**1/2)') 
ALPHA2 = ALPHA * 1055.D0 
WRITE (IWTTY,220) ALPHA,ALPHA2 
220 FORMAT (' FIRE GROWTH CONSTANT = ',F12.4,' BTU/(SEC**3)'/ 
& : (',E12.4,' JOULES/SEC**3)') 
WRITE (IWTTY, 400) 
READ (IRTTY,410) M 
WRITE (IWTTY, 230) 
230 FORMAT (/'O FOR TEMPERATURE ACTUATED DETECTOR: ' ) 
TACT = TACT - 460.D0 
TACT2 = (5.D0/9.D0)*(TACT-32.D0) 
WRITE (IWTTY,240) TACT, TACT2 
240 FORMAT ('0 ACTIVATION TEMPERATURE = ',F6.1,' DEGREES F (' ,F6.1, 
& ' DEGREES C)') 
T=T / 60.D0 
WRITE (IWTTY,250) T 
250 FORMAT ('0 TIME TO ACTIVATION = ',F8.2,' MINUTES' ) 
QD2 = QD * 1.055D0 
WRITE (IWTTY,260) QD,QD2 
260 FORMAT ('0 HEAT RELEASE RATE = ',E12.4,' BTU/SEC'/ 
& ; (',E12.4,' KILOJOULES/SEC)') 
WRITE (IWTTY,400) 
READ (IRTTY,410) M 
WRITE (IWTTY, 270) 
270 FORMAT (//'O FOR RATE OF RISE ACTUATED DETECTOR: ') 
ROR2 = ROR * (5.D0/9.DO) 
WRITE (IWTTY,290) ROR,ROR2 
290 FORMAT ('0 ACTIVATION RATE OF RISE = ',F6.2,' DEGREES F/MIN (', 
& F6.2,' DEGREES C/MIN)') 
TR = TR / 60.D0 
WRITE (IWTTY,250) TR 
QDR2 = QDR * 1.055D0 
WRITE (IWTTY,260) QDR,QDR2 
WRITE (IWTTY, 400) 
READ (IRTTY,410) M 


ELSE 

C 

C METRIC (ENGLISH) OUTPUT 

G 
HF2 = HF * (1.D0/0.3048D0) 
ZF2 = ZF * (1.D0/0.3048D0) 


WRITE (IWTTY,300) HF,HF2,ZF,ZF2 

300 FORMAT ('0O CEILING HEIGHT = ',F6.2,' METERS (',F6.2,' FEET)'/ 
& 5X, 'DETECTOR SPACING = ',F6.2,' METERS (',F6.2,' FEET)') 
RTI2 = RTI * ((1.D0/0.3048D0)**0. 5D0) 
WRITE (IWTTY,320) RTI,RTI2 

320 FORMAT ('0 DETECTOR RTI = ',F6.1,' (M-SEC)**1/2 (',F6.1,' (FT-S 
&EC)**1/2)') 
ALPHA2 = ALPHA * (1.D0/1055.D0) 
WRITE (IWTTY,330) ALPHA,ALPHA2 
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330 


340 


350 


360 


370 


C 


C 


FORMAT ('0 FIRE GROWTH CONSTANT = ',E12.4,' JOULES/(SEC**3)!/ 
& ' (',E12.4,' BTU/SEC**3)') 
WRITE (IWTTY, 400) 

READ (IRTTY,410) M 

WRITE (IWTTY, 230) 

TAGE = TACT —' 273). D0 

TACT2 = (9.D0/5.D0) * TACT + 32.D0 

WRITE (IWTTY,340) TACT, TACT2 

FORMAT ('0 ACTIVATION TEMPERATURE = ',F6.1,' DEGREES C (',F6.1, 
& ' DEGREES F)') 

T=T / 60.D0 

WRITE (IWTTY, 350) T 

FORMAT ('0 TIME TO ACTIVATION = ',F8.2,' MINUTES’ ) 

QD = QD / 1000.D0 

QD2 = QD * (1.D0/1.055D0) 

WRITE (IWTTY,360) QD,QD2 

FORMAT ('0 HEAT RELEASE RATE = ',E12.4,' KILOJOULES/SEC' / 
& : (',E12.4,' BTU/SEC)') 

WRITE (IWTTY, 400) 

READ (IRTTY,410) M 

WRITE (IWTTY, 270) 

ROR2 = ROR * (9.D0/5.D0) 

WRITE (IWTTY,370) ROR,ROR2 

FORMAT ('0 ACTIVATION RATE OF RISE = ',F6.2,' DEGREES C/MIN (', 
& F6.2,' DEGREES F/MIN)') 

TR = TR / 60.D0 

WRITE (IWTTY,350) TR 

QDR = QDR / 1000.D0 

QDR2 = QDR * (1.D0/1.055D0) 

WRITE (IWTTY,360) QDR,QDR2 

WRITE (IWTTY, 400) 

READ (IRTTY,410) M 

ENDIF 

WRITE (IWTTY, 420) 

READ (IRTTY,410) M 

iP C(MeEO OY!) ORM CMaEO. "y')))) COTO! 5 

STOP 'PROGRAM COMPLETED! 

CONTINUE 

WRITE (IWTTY, 390) 

FORMAT (' **%:%2%%> ERROR IN DETACT ROUTINE <*%%%%%%! ) 

FORMAT (//' <RETURN> TO CONTINUE! ) 

FORMAT (Al) 

WRITE (IWTTY, 420) 

FORMAT (' TRY ANOTHER CASE (Y/N) ?'/) 

READ (IRTTY,410) M 

IF ((M.EQ.'Y').OR.(M.EQ.'y')) GO TO 5 

STOP 'PROGRAM ABORTED' 

END 


SUBROUTINE NWIN (C,P,IERR, IWTTY) 


C THIS IS ROUTINE NWIN. IT SOLVES THE FIXED TEMPERATURE DETECTOR 
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C ACTIVATION TIME EQUATION USING A NEWTON-RAPHSON METHOD. 
G 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 
PO = 0.1D0 
TOL = 0.00001D0 
NO = 1000 
te il 
TERR = 0 
10 CONTINUE 
IF (I.LE.NO) THEN 
IF (P0.GT.50.D0) THEN 
X = 0.0D0 
ELSE 
X = DEXP(-PO) 
ENDIF 
FX = PO+X-C 
FPMX = 1.0D0 - X 
IF (FPMX.LT.0.00000001D0) GO TO 30 
P = PO - (FX/FPMX) 
IF (ABS(P-PO).LT.TOL) THEN 
TERR = 0 
RETURN 
ELSE 
IS ie cb i 
PO = P 
ENDIF 
GO TO 10 
ENDIF 
TERR = 1 
WRITE (IWTTY,20) I 
20 FORMAT (' NEWTON-RAPHSON FAILED AFTER NO ITERATIONS, NO = ',I4) 
RETURN 
30 WRITE (IWTTY,40) 
40 FORMAT (' SLOPE OF EQUATION TOO CLOSE TO ZERO FOR '/ 


& ' NEWTON-RAPHSON METHOD. '/ 
& ' ERROR RETURN' ) 
TERR = 1 
RETURN 
END 
C 
SUBROUTINE BISECT (D1,D2,ROR,P,IERR, IWTTY) 
@ 
C THIS IS ROUTINE BISECT. IT SOLVES THE RATE OF RISE DETECTOR ACTIVATION 
C TIME EQUATION USING A BISECTION METHOD. 
C 
IMPLICIT DOUBLE PRECISION (A-H,0-Z) 
TERR = 0 
TOL = 0.00001D0 
NO = 1000 
A = 0.0D0 
B = 1000.0D0 
RLMT = TOL / 20.0 


40 


10 


20 


30 


CONTINUE 

IF ((D2*B).GT.50.D0) THEN 
X = 0.0D0 

ELSE 

X = DEXP (-D2*B) 

ENDIF 

FXB = D1*B**0.25D0-D1*B**0.25*X-ROR 
IF (FXB.LT.0.0D0) THEN 
A=B 
B = B + 500.D0 

GO TO 10 

ENDIF 

= il 

CONTINUE 

IF (I.LE.NO) THEN 

P = A+(B-A)/2.0D0 

IF ((D2*P).GT.50.D0) THEN 

X = 0.0D0 

ELSE 

X = DEXP (-D2*P) 

ENDIF 

FX = D1*P**0.25D0-D1*P**0.25D0*X-ROR 

IF (((FX.GT.-RLMT).AND. (FX.LT.RLMT)).OR. 
&  (((B-A)/2.D0).LT.TOL)) THEN 

TERR = 0 

RETURN 

ELSE 

f= P+ 

IF ((D2*A).GT.50.D0) THEN 

X = 0.0D0 

ELSE 

X = DEXP (-D2*A) 

ENDIF 

FXA = D1*A**0.25D0-D1*A**0.25D0*X-ROR 

IF ((FXA*FX).GT.0.0D0) THEN 

A=pP 

ELSE 

B=pP 

ENDIF 

ENDIF 

GO TO 20 

ENDIF 

IERR = 1 

WRITE (IWTTY,30) I 

FORMAT (' BISECTION ROUTINE FAILED AFTER NO ITERATIONS, NO = ',I4) 
RETURN 

END 
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APPENDIX B. DETACT-T2 
BASIC VERSION 
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10 PRINTS DETACT- 12)  WVERSTON G2 0u: 

20 PRINT “WRITTEN BY D.W. STROUP (1985)" 

SOP RUN Tsai 

40 PRINT "CONTRIBUTION OF THE" 

50 PRINT "NATIONAL BUREAU OF STANDARDS (U.S.)." 
60 PRINT "NOT SUBJECT TO COPYRIGHT. " 

OME RAUNG ie 

80 PRINT "CALCULATES DETECTOR ACTUATION TIME" 
90 PRINT "BELOW UNCONFINED CEILINGS WITH " 

100 PRINT "TIME SQUARED FIRE GROWTH RATES " 


110) BRON" 

120 REM 

130 REM INPUT: 

140 REM J - UNITS CODE (1 OR 2) 

150 REM 1 - INPUT DATA IN ENGLISH UNITS 
160 REM 2 - INPUT DATA IN METRIC UNITS 
170 REM TAMB - AMBIENT TEMPERATURE 

180 REM RTI  - DETECTOR RESPONSE TIME INDEX 

190 REM TACT - DETECTOR ACTIVATION TEMPERATURE 
200 REM ROR’ - DETECTOR ACTIVATION RATE OF RISE 
210 REM HF - CEILING HEIGHT 

220 REM ZF - DETECTOR SPACING 

230 REM M - GROWTH RATE FACTOR CODE, CHARACTER (S, M, F, U, OR 0) 


240 REM ALPHA - FIRE GROWTH RATE CONSTANT (FOR M = 0) 
250 REM 
260 REM OUTPUT: 


ZO REM ies - TIME OF ACTIVATION FOR FIXED TEMPERATURE DETECTOR 
280 REM QD - HEAT RELEASE RATE AT TIME OF ACTIVATION, T 
290 REM TR - TIME OF ACTIVATION FOR RATE OF RISE DETECTOR 
300 REM QDR  - HEAT RELEASE RATE AT TIME OF ACTIVATION, TR 
310 REM IERR - ERROR CODE (0 OR 1) 

320 REM 0 - SUCCESSFUL 

330 REM 1 - UNSUCCESSFUL 

340 REM 

350 CLEAR 

360 REM 

370 REM ENTER INPUT DATA 

380 REM 


390 READ GE#,CPE##, RHOE# 

400 DATA 32. 2i#,0.24#,0.0735# 

410 READ GM#f,CPM#, RHOM# 

420 DATA 9.8#,1.0035#,1.1768# 

430 IERRZ = 0 

440 PRINT " ENTER: 1 FOR ENGLISH UNIT INPUT" 
450 PRINT " 2 FOR METRIC UNIT INPUT" 
460 INPUT "SELECTION = ",JZ 

470 IF JZ<>1 GOTO 740 

480 REM 

490 REM ENTER DATA IN ENGLISH UNITS 

500 REM 
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510 INPUT "AMBIENT TEMPERATURE <DEGREES F> = '",TAMB# 

520 INPUT "DETECTOR RESPONSE TIME INDEX, RTI <(FT-SEC) (1/2)> = '",RTI# 
530 IF RTI#<.0000001# THEN LET RTI# = .0000001# 

540 RTIR# = RTI# 

550 INPUT "DETECTOR ACTIVATION TEMPERATURE <DEGREES F> = '', TACT 

560 INPUT "DETECTOR RATE OF TEMPERATURE RISE <DEGREES F/MIN> = ',ROR# 
570 INPUT ''ROOM CEILING HEIGHT <FEET> = ",HF# 

580 INPUT "DETECTOR SPACING <FEET> = "', ZF# 

590 PRINT "ENTER: S FOR SLOW FIRE GROWTH RATE," 

600 PRINT "' M FOR MEDIUM FIRE GROWTH RATE," 

610 PRINT '' F FOR FAST FIRE GROWTH RATE," 

620 PRINT '' U FOR ULTRAFAST FIRE GROWTH RATE, OR'' 

630 PRINT '' O FOR OTHER FIRE GROWTH RATE" 

640 INPUT ''SELECTION = ",M$ 

650 IF (M$="'s'" OR M$="S'') THEN LET ALPHA# = .00277778# 

660 IF (M$="'m'' OR M$="M'"') THEN LET ALPHA# = .01111111# 

670 IF (M$="£" OR M$="F'') THEN LET ALPHA# = .04444445} 

680 IF (M$="'u'' OR M$="U") THEN LET ALPHA# = .17777778# 

690 IF (M$="0'"' OR M$="'0'"') THEN INPUT "FIRE GROWTH RATE <BTU/SEC 3> = '', ALPHA# 
700 GOTO 940 

710 REM 

720 REM ENTER DATA IN METRIC UNITS 

730 REM 

740 INPUT "AMBIENT TEMPERATURE <DEGREES C> = ',TAMB# 

750 INPUT "DETECTOR RESPONSE TIME INDEX, RTI <(M-SEC) (1/2)> = ",RTI# 
760 IF RTI#<.00000014 THEN LET RTI# = .0000001# 

770 RTIR# = RTI# 

780 INPUT ''DETECTOR ACTIVATION TEMPERATURE <DEGREES C> = ", TACT# 

790 INPUT ''DETECTOR RATE OF TEMPERATURE RISE <DEGREES C/MIN> = "', ROR# 
800 INPUT 'ROOM CEILING HEIGHT <METERS> = ",HF# 

810 INPUT ''DETECTOR SPACING <METERS> = '',ZF# 

820 PRINT "ENTER: S FOR SLOW FIRE GROWTH RATE," 

830 PRINT '' M FOR MEDIUM FIRE GROWTH RATE," 

840 PRINT " F FOR FAST FIRE GROWTH RATE," 

850 PRINT '' U FOR ULTRAFAST FIRE GROWTH RATE, OR'' 

860 PRINT '"' O FOR OTHER FIRE GROWTH RATE" 

870 INPUT "SELECTION = '"',M$ 

880 IF (M$="'s'' OR M$="S") THEN LET ALPHA# = 2.9305556# 

890 IF (M$="'m'' OR M$="M'') THEN LET ALPHA# = 11.722222# 

900 IF (M$=""£" OR M$="F'') THEN LET ALPHA# = 46.888889# 

910 IF (M$="u"' OR M$="U'"') THEN LET ALPHA# = 187.55556# 

920 IF (M$="0"' OR M$="0'') THEN INPUT "FIRE GROWTH RATE <J/SEC 3> = '",ALPHA# 
930 REM 

940 REM Cars aes seade Sets aes tetas dead CALCULATIONS PT aoe do dd don dn ca nahn do gb ag do dint ohn ale at niet aieurabow se oe salpa oa des des ds dp 
950 REM 

960 Ri# = .5S# * SQR(2#) * ZF# 

970 ROH# = R# /HF# 

980 ROR# = ROR# / 60# 

990 IF JZ<>1 GOTO 1050 

1000 TAMB# = TAMB# + 460} 


1010 TACT# 


TACT# + 460# 


1020 Ad = GE##/ (CPE#*TAMBi#*RHOE# ) 
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GH = GE# 
GOTO 1090 
TAMB# = TAMBi# + 273+ 
TACT# = TACT# + 273# 


Aft = GM##/ (CPM#*TAMB##*RHOM#* 1 000#) 


GH = GM# 


TOTS2#=A# (-14/54#)*ALPHA# (-14/5i)*HFH (44/54) 
DLTODL#=A# (24#/5i#)*(TAMBH#/Gi)*ALPHAH (2d/5#)*HF# (-34/ 54) 
UOUS2#=A#” (14/54) *ALPHA# (14/54) *HFH (1H / 54) 


DELTD# = TACT# - TAMB# 


IF (ROH#>.3) THEN UODLTH#=.59*ROH# (-.63#) ELSE UODLTH#=3.87#/(SQR(9.115#)) 


TS2F# = .954f * (1i# + ROH#) 
P# = TOTS2#*(. 188+. 3134*ROH# ) 


A2#=(4#/ 34#)*DLTODL#*UOUS2# (- . 54) *UODLTH# (- . 54) *RTI#/ Pit 


CH = 14 + DELTD# / A2# 
GOSUB 2650 
IF IERR%Z=1 GOTO 2550 


DELTS2#=(4#/34#)*UOUS2# (-. 54) *UODLTH# (- . 54) *RTI#* Yi / Pit 
TS2i#=. 954% ( 1#+ROH#)+(. 1884+. 3134*ROH#)*DELTS2#° (34 / 4H) 


TS2# = TS2# + TS2Fi# 


TH=TS2#*AH (- 14/54) *ALPHAH (- 14/54) HEH (4d / 544) 


QD# = ALPHA * T/°2! 


IF RTIR#<=2/ THEN LET RTIR# = 2# 


Di ft=( 44/34) *DLTODL#/ Pi 


D2#=( 34/44) *UOUS2# . 54 UODLTH# . 54#*( 141 /RTIR#) * PH 


GOSUB 2950 
IF IERRZ=1 GOTO 2550 


TS2Rit=. 954d (14+ ROH#)+(. 1884+. 3134#*ROHH#)*DELTS2# (34/44) 


TS2R# = TS2R# + TS2F# 


TR#=TS2RH*AH (-14/54)*ALPHAH (-14/54)*HF# (4i/ 54) 


QDR# = ALPHA# * TR# 2! 
ROR# = ROR# * 60# 

REM 

REM PRINT OUT RESULTS 
REM 

PRINT "RESULTS:"' 


IF J%<>1 GOTO 1960 
REM 
REM ENGLISH (METRIC) OUTPUT 
REM 


HF2# = HF# * .3048 

ZF2# = ZF# * .3048 

PRINT " "' 

PRINT '' CEILING HEIGHT = '";CSNG(HF#);'' FEET" 

PRINT " ("'; CSNG(HF2#) ;"' METERS)" 
PRINT "' "' 

PRINT '"' DETECTOR SPACING = ";CSNG(ZF#);"' FEET" 
PRINT " ('";CSNG(ZF2#);'" METERS)" 
PRINT" "! 

RTI2# = RTI# * (.3048# .5#) 














PRINT " RTI = ";CSNG(RTI#);'' (FT-SEC)~0.5";" (''s;CSNG(RTI2#);" (M-SEC) 0.5) 


PRINT " 
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ALPHA2# = ALPHA# * 1055# 

PRINT " FIRE GROWTH RATE = '";CSNG(ALPHA#);'' BTU/(SEC 3)" 

PRINT " (''; CSNG(ALPHA2#);'' J/(SEC 3))" 

PRINT "' "' 

INPUT ''<RETURN> TO CONTINUE "';0$ 

PRINT " " 

PRINT " "! 

PRINT '' FOR TEMPERATURE ACTUATED DETECTOR:"' 

TACT# = TACT# - 460+ 

TACT2# = (5#/94) * (TACTH# - 32#) 

PRINT " "' 

PRINT " ACTIVATION TEMPERATURE = ";CSNG(TACT#);'' DEGREES F" 
PRINT " ('';CSNG(TACT2#) ;'' DEGREES C)'' 
PRINT " "' 

TH = TH / 604 

PRINT '' TIME TO ACTIVATION = ";CSNG(T#);'' MINUTES" 

PRINT " '"' 

QD2# = QD# * 1.055# 

PRINT " HEAT RELEASE RATE = ";CSNG(QD#);'' BTU/SEC" 

PRINT " (";CSNG(QD2#);" kJ/SEC )" 

PRINT " "' 

INPUT '"<RETURN> TO CONTINUE ";0$ 

PRINT " "' 

PRINT " "! 

PRINT "' FOR RATE OF RISE ACTUATED DETECTOR:" 

ROR2# = ROR# * (5i#/9#) 

PRINT " "' 

PRINT '' ACTIVATION RATE OF RISE = ";CSNG(ROR#);'! DEGREES F/MIN" 
PRINT " (''; CSNG(ROR2#);'' DEGREES C/MIN)" 
PRINT "' "' 

TR = TRi# /60# 

PRINT '' TIME TO ACTIVATION = '";CSNG(TR#);'' MINUTES" 

QDR2# = QDR# * 1.0554 

PRINT " "' 

PRINT '' HEAT RELEASE RATE = ''sCSNG(QDR#);'' BTU/SEC" 

PRINT "' (";CSNG(QDR2#);" kJ/SEC )" 

PRINT " "' 

PRINT " "' 

INPUT '' TRY ANOTHER RUN (YES/NO) ';0$ 

IF (O$="Y" OR O$="'y'"') GOTO 120 

GOTO 2530 

HF2# = HF * (1!/.3048) 

ZF2# = ZF# * (1!/.3048) 

REM 

REM METRIC (ENGLISH) OUTPUT 

REM 
PRINT 
PRINT '' CEILING HEIGHT = '";CSNG(HF#);'' METERS" 
PRINT " ("3 CSNG(HF2#);"' FEET)" 
PRINT "' " 

PRINT '' DETECTOR SPACING = '';CSNG(ZF#);'' METERS" 
PRINT "' (";CSNG(ZF2#);'' FEET)" 


- 
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PRINT 
RTI2# 
PRINT 
PRINT 


ALPHA2}} 


PRINT 
PRINT 
PRINT 
INPUT 
PRINT 
PRINT 
PRINT 
PRINT 
TACT# 


RTI# * ((11/.3048)~ . 5#) 


RTI = '"sCSNG(RTI#);" (M-SEC)~0.5"'s" (''sCSNG(RTI2#);"" (FT-SEC)~0.5) 


= ALPHA# * (11/1055!) 
FIRE GROWTH RATE = '";CSNG(ALPHA#H);" J/(SEC 3)" 
("'; CSNG(ALPHA2#) ;" BTU/(SEC 3))" 


"<RETURN> TO CONTINUE '"3;0$ 
Ue 


FOR TEMPERATURE ACTUATED DETECTOR:" 
it 


TACT# - 273# 


TACT2# = (Of/5i#) * TACTH# + 32+ 


PRINT 
PRINT 


"0 


ACTIVATION TEMPERATURE = '";CSNG(TACT#);'' DEGREES C" 
("'; CSNG(TACT2#);" DEGREES F)" 


PRINT 
TH = TH / 60# 


PRINT 


TIME TO ACTIVATION = ";CSNG(T#);'' MINUTES" 


PRINT "' "' 
QD# = QDi# / 1000} 
QD2# = QD# * (1!/1.055) 


PRINT 
PRINT 


HEAT RELEASE RATE = '';CSNG(QD#);" kJ/SEC" 
(''s CSNG(QD2#);'' BITU/SEC )" 


PRINT 
INPUT ''<RETURN> TO CONTINUE"; O$ 


PRINT 
PRINT 
PRINT 
ROR2#} 
PRINT 
PRINT 
PRINT 


at 


FOR RATE OF RISE ACTUATED DETECTOR:" 
ROR# * (94/5) 
0 


ACTIVATION RATE OF RISE = ';CSNG(ROR#);'' DEGREES C/MIN" 
(""; CSNG(ROR2#);'' DEGREES F/MIN)" 


BRANT 
TR# /60# 


TR# = 
PRINT 
PRINT 


TIME TO ACTIVATION = ";CSNG(TR#);'' MINUTES" 


QDR# = QDR# / 1000# 


QDR2i! 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


QDR# * (11/1.055) 
HEAT RELEASE RATE = '";CSNG(QDR#);" kJ/SEC" 
('";CSNG(QDR2#);"" BTU/SEC )" 


"9 
Tt] 


TRY ANOTHER RUN (YES/NO) '';0$ 


IF (O$="Y" OR O$="y'"') GOTO 120 


END 
REM 


REM xvekoevekseksex% ERROR HANDLING ROUTINE *okvecvesesesedeses 


REM 


PRING. 
PRINT " !!!!! ERROR IN DETACT-T2 ROUTINE I!!!1!" 
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PRINT " !!!!! PROGRAM RUN ABORTED eae 
PRIND 1" 

INPUT '' WOULD YOU LIKE TO TRY AGAIN (YES/NO) ?'';0$ 
IF (O$="Y" OR O$="y") GOTO 120 


END 

REM 

REM ** NEWTON RAPHSON ROUTINE TO FIND ACTIVATION ** 
REM ** TIME OF FIXED TEMPERATURE DETECTOR. wx 
REM 

PO# = .1# 

TOL# = .00001# 

NOZ = 1000 

IZ =1 

IERRZ = 0 

REM 

IF IZ>NOZ GOTO 2870 


IF (PO#>50#) THEN LET X#=O/f ELSE X#=EXP(-PO#) 
FX#=PO/t + Xi## - CH 

FPMX# = 1# - Xi 

IF (FPMX#<.00000001#) GOTO 2900 
P2# = PO# - (FXi#/FPMX#) 

IF (ABS(PO#-P2#)>TOL#) GOTO 2840 
ITERR% = 0 

Yi = P2+ 

RETURN 

1S Wy 45 1 

PO# = P2# 

GOTO 2730 

IERRZ = 1 


PRINT "NEWTON RAPHSON FAILED AFTER NO ITERATIONS, NO = ",1Z% 


RETURN 
PRINT "SLOPE OF EQUATION TOO CLOSE TO ZERO!!!" 
PRINT “ERROR IN NEWTON RAPHSON ROUTINE" 


IERRZ = 1 

RETURN 

REM 

REM ** BISECTION ROUTINE TO FIND ACTIVATION ** 
REM ** TIME OF RATE OF RISE DETECTOR. we 
REM 

TERRZ = 0 

TOL# = .00001# 

NOZ = 1000 

AA = Of 

Bi = 1000# 


RLMT# = TOL# / 20# 
REM 


IF (D2#*Bi#)>50# THEN Xi# = Of ELSE X# = EXP(-D2/*B#f) 
FXB# = Dl#t * Bi ~ .25# - Diff * Bi ~ .254# * X# - ROR# 
IF FXBi#>=0# GOTO 3110 

AAi# = Bit 

Bi = Bit + 500# 

GOTO 3040 
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1 = il 
REM 

IF I>NO GOTO 3260 

P3# = AA + (Bi - AAH) / 24 

IF (D2#*P3#)>50# THEN X# = Off ELSE Xi## = EXP (-D2#*P3#) 
FX# = Dl# * P3# ~ .25# - D1# * P3H# ~ .254 * Xt - ROR# 
DELTS2# = P3# 

IF ((FX#>(-RLMT#)) AND (FX#<RLMT#)) THEN RETURN 

IF ((BiH-AA#) /2#)<TOL# THEN RETURN 


REM 
1 eB Ib4 eb il 
IF (D2#*AA#)>504 THEN X# = Off ELSE X# = EXP(-D2/*AA#) 


FXA# = D1i# * AA# ~ .25# - Didt * AA ~ .254 * Xf - ROR# 

IF (FXA#*FX#)>0/ THEN AA/# = P34 ELSE BH = P3# 

GOTO 3120 

TERRZ = 1 

PRINT "BISECTION ALGORITHM FAILED TO CONVERGE AFTER ",1Z,'' ITERATIONS." 
RETURN 
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APPENDIX C. 


DETECTOR ACTIVATION TIME TABLES: 
FIXED TEMPERATURE 
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